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Abstract. NGC 4258 is well known for its 22 GHz water maser features that 
appear to reside in a thin Keplerian disk orbiting at i?o = 0.13 pc. The influence 



> , 

of gravitational redshift on the level of maser saturation and the line strength 

and width of systemic water maser emission is studied. The effective gravitational 

velocity shift along the gain path through the disk amounts to about twice the 

l/^ , thermal velocity width of H2O molecules at 500 K for a masing ring of thickness 

, 0.27?o- This GR induced frequency shift has a sign opposite to the shift caused 

by the projected velocity difference vector of Keplerian motion on the blue side 

q . of the systemic maser emission and has the same sign on the red systemic side. 

S_i ■ 

+^ ■ Consequently, the negative optical depth that a line photon sees as it traverses 

C/3 . 

the masing disk along a path of minimal Keplerian velocity gradient is larger on 

the blue side. Note that one is not concerned here with an overall shift in the 
• i-H , 

maser spectrum caused by the fact that the masing ring as a whole is located in 
the black hole's potential well, i.e., it is the change in the gravitational potential 
along the gain path that matters. 

It is found, for a Schwarzschild metric, that saturation is enhanced (sup- 
pressed) for systemic blue-shifted (red-shifted) maser features by gravitational 
redshift of line photons. This asymmetry should be observable for NGC 4258. 
That is, approaching peak flux systemic features are, on average, i.e., over a pe- 
riod of several years, brighter, due to saturation, and wider (~ 2 km/s), due to 
rebroadening, than their receding counterparts (~ 1 km/s). The effects are con- 
sistent with observations that span several years, i.e., over periods longer than 
needed for disk material to move through the systemic emission region. A com- 
pletely definitive confirmation of this effect can follow by combining all available 
systemic 22 GHz data over the past two decades. 

keywords, galaxies: individual: NGC 4258 - ISM: molecules - molecular processes 
- radiative transfer 
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1. Introduction 

Bright 22 GHz water masers have been detected toward NGC 4258 (M106) and other 
galaxies (Nakai et al. 1993). For the AGN NGC 4258 the detection of both high- velocity, 
~ —400 and ~ 1400 km/s, and systemic, ~ 500 km/s, emission (LSR) have lead to the 
picture of a (warped by 6°) Keplerian accretion disk that rotates, at a distance of about 
0.1 pc, with <~ 900 km/s around a massive (M ss 3.5 x 10 7 M Q ) black hole (Moran et al. 
1995). The overall appearance of the spectra is spiky and individual maser features have 
FWHM widths of 1 - 2 km/s (Grcenhill ct al. 1995, 1996). The excitation (pumping) of 
the masers involves the the efficient formation of H 2 at temperatures in excess of a few 
hundred Kelvin at H 2 densities n > 10 8 cm~ 3 (Neufeld et al. 1994). The appearance of 
high-velocity and systemic maser emission, as well as the presence of absorption in the 
central feature, has been explained by Watson & Wallin (1994) and Wallin ct al. (1998) 
for unsaturated amplification. Here, the effects of saturation and gravitational redshift 
caused by the central black hole are studied. 

2. Model Description 

A thin Keplerian accretion disk, viewed edge-on, with a height H = 10 15 cm and radial 
coordinate r, is adopted as in Wallin et al. (1998). For the impact parameter b one can 
write the optical depth of unsaturated masing gas as 



where a is the (constant) maser opacity, the frequency obeys v = isq(1 — 5v/c) for the 
rest frequency vq and the equivalent velocity difference 8v. Gravitational redshift adds a 
contribution Svqr = cr g (l/r — 1/r') for the Schwarzschild radius r g — 2M in geometric 
units and the black hole mass M. Keplerian motion follows 6(GM) 5 /r 1,5 . Radial gas 
motion (inward drift < 0.007 km/s) is ignored in this work (Moran et al. 1995). One finds 
r g = 1.1 x 10 cm. The Doppler width vu = (v^ h + U( U ) ' 5 contains the thermal width 
v t h = (2kT/m) ' 5 = 0.7 km/s of the H2O gas molecules with mass m at T = 500 K as well 
as turbulent motions Vt u . It is assumed that the systemic masering water molecules are 
located at Ro = 0.13 pc over a length scale Ar = 0.2i?o, while the entire masing disk runs 
to 2i? (Greenhill et al. 1995; Wallin ct al. 1998, 1999). The Keplerian rotational velocity 
ranges between v rot ~ 800 — 1100 km/s from the outside to the inside of the masing disk 
(Hashick ct al. 1994). The unsaturated gain 7 at line center is taken as ln{^) ~ 11 at 
LSR velocities of ±20 km/s (with respect to systemic), which yields ln(-f) ~ 24 at about 
±900 km/s. These numbers roughly reproduce a typical observed flux density of about 2 
Jy with a temperature for the (unresolved) background continuum source of T c ~ 10 7 K 
for b ~ (Wallin et al. 1998) and T c ~ 10 K for b > 1 (almost self-amplification; Watson 
& Wallin 1994), where b is normalized to Ro- It is assumed that a layer of non- inverted 
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H2O molecules outside or inside of the masing gas, but with the same Doppler width 
and a large (> 25 at b — 0) optical depth, is reponsible for the absorption around the 
systemic velocity (Wallin et al. 1998). The distance to NGC 4258 is 6.4 Mpc (Miyoshi et 
al. 1995). The measured systemic velocity of 472 ± 4 km/s derives from Ha observations 
(Cecil et al. 1992). The results presented below do not depend strongly on the assumed 
model parameters. 

3. Asymmetric Systemic Emission: Gravitational Redshift and Saturation 

The rate of de-excitation due to stimulated emission, R s , is, for an arbitrary point q, 

R s (q) = J dQR s (q,n) = J dvc 2 : A/8nhu% J (j){v)I v (£l)dQ. (2) 

and depends on the magnitude of the velocity gradient along the line of sight since this 
determines which path segment contributes to the frequency and solid angle integrals. In 
this, A is the Einstein A coefficient of the 22 GHz maser transition, Vi is the solid angle 
and <f> is the absorption profile centered at vq. 

Saturation for the 22 GHz maser line occurs when stimulated emission competes with 
the effective decay rate, T ~ 1 s _1 (~ 0.1 s _1 only if infrared transitions are very optically 
thick in all directions), of masing states due to all other processes. The rate for velocity 
relaxation, R v (set by the lifetimes of infrared decays), is at least as rapid and R v ~ 2s _1 
at n = 10 9 cm -3 and 500 K (Anderson & Watson 1993). Rebroadening of maser lines is 
important when R s exceeds R v (Nedoluha & Watson 1991). For NGC 4258, observations 
(including Rq and Ar) indicate that R s ~ T (Wallin et al. 1998) and, hence, a significant 
increase in R s would strongly impact the saturation and line shapes of maser features. 

One obtains for the step function 9 and vu = vqVd / c that 

<f>(v) » [6(vo - v D ) - 6{v + v D )\/2v D (3) 

and, to a good approximation, 

R s (q) = J dflR s (q - L[g], fi) e _a dp (4) 

for the velocity gradient g = dv/dp along a path parameterized through p(q) and L = 
Vo/\g\ the local Sobolev length at line center. If one dominant beam exists, the problem 
becomes effectively one-dimensional, and one can concentrate on a single direction and 
work only with R s (q, O). When one applies the variational principle to find the change 
in the Sobolev length that results from Kcplerian motion and gravitational redshift, the 
result is 

SR S (Q) w R s (Q)e TKSL/L (5) 

with SL = L[KGR] — L[K], in which K indicates Kcplerian motion and KGR combined 
Keplerian motion and gravitational redshift. Evaluation of 5L is straightforward through 

Svqr = cr g (l/r - 1/r') km/s, (6) 



4 



Marco Spaans: Water Masers in NGC 4258 - GR and Saturation 
Sv K = (GMf ^bRoil/r 15 - 1/r' 1 - 5 ) km/s, 



(7) 



and indicates that gravitational redshift compensates the Doppler shift that results from 
projected Keplerian motion on the approaching side of the accretion disk. The converse 
is true for the receding side. Hence, for the systemic maser features saturation is both 
enhanced and suppressed. It can be concluded immediately that blue-shifted systemic 
maser features should be, on average, brighter than their red-shifted counterparts. The 
effect of gravitational redshift is quite small for the high-velocity features since only a 
narrow range in \r — r'\ is probed. Also, special relativistic effects are neglected since 
changes in the 7 factor are extremely small due to the required velocity coherence over 
a line width of 1-2 km/s. 

Over a scale of Ar at b = 0, gravitational redshift amounts to Svqr = 1-4 km/s. It is 
further assumed that v tu — v s with the sound speed v s — 1.3 km/s at 500 K. A Monte 
Carlo method, see Spaans & van Langevelde (1992) for details, is used to evaluate 5R S . 
It is found for the adopted parameters that the average Sobolev length L increases for 
the approaching systemic features by about 30%, over Ar. This causes the decay rate 
due to stimulated emission to increase from R s ~ 1 s _1 to R s {—) ~ 9 s _1 for blue-shifted 
maser features over the observed LSR velocity range, and with respect to the originally 
adopted gain 7. Conversely, one still finds R s {+) ~ 0.8 s _1 for red-shifted maser emission 
along the longest path, i.e., the path that originates at about a velocity width from 
line center. Hence, since the adopted parameters for Rq, Ar, M and 7 are motivated 
by observations, it can be concluded generally that blue systemic features should be 
impacted more strongly by saturation than red systemic features. Note in this that a 
saturated maser source becomes a "linear convertor" for R s (— ) >T/h with h <~ 1/2 and 
thus that the brightest blue systemic features should have peak flux densities that are 
larger than their red counterparts by a factor of F ~ T / (hR s (+)) ~ 2.5, consistent with 
observations of peak flux densities (~ 5 Jy) of individual blue systemic maser features. 

For the blue-shifted systemic maser features in NGC 4258 one expects, in the (almost) 
saturated limit, a typical FWHM line width of Av < 1.67vd = 2.5 km/s, but in any case, 
at least, Av > v D = 1.5 km/s (Emmering & Watson 1994; Nedoluha & Watson 1991), 
consistent with observations (Moran et al. 1995; Greenhill et al. 1995). The FWHM width 
of a typical red-shifted systemic maser feature decreases only modestly with R s and is 
~ 0.6vd = 1.0 km/s. Hence, the brightest red-shifted systemic features should be, on 
average, narrower by a factor of ~ 2 than the brightest blue-shifted systemic features. 

It is concluded that turbulence, which is likely present in accretion disks (Afshordi 
et al. 2004), need not be very large, i.e., be nearly subsonic, in order to explain the 
typical FWHM widths, 1-2 km/s, of emission spikes in observed spectra of NGC 4258. 
Supersonic turbulence decays rapidly, while subsonic turbulence is easier to maintain. 
This is consistent with observations that the disk is not much thicker than <~ 10 15 cm 
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(Moran et al. 1995) given that turbulence alone can support a disk with height H ~ 
VtuRo/ v rot — 5 x 10 14 cm. Finally, this result differs somewhat from that of Wallin et al. 
(1998), who argue for a larger turbulent velocity in an unsaturated medium to reproduce 
the observed maser feature FWHM widths, but agrees with the tilted disk model of 
Wallin et al. (1999). 

4. Observational Comparison and Caveats 

Observations indicate that blue-shifted systemic emission features tend to be brighter 
and wider, on average, than receding systemic features (Hashick et al. 1994, their figure 
2; Wallin et al. 1998, their figure 3). It is stressed here that the data displayed in figure 
2 of Hashick et al. (1994) is cumulative over a period of 5 years, 1987-1991. Over such a 
period any fluctuation in the ambient conditions caused, e.g., by turbulence, would travel 
from about -25 to +25 km/s in the systemic emission region, approximately 0.005 pc, if 
the fluctuation moves with the bulk Keplerian motion. That is, such a fluctuation would 
be amplified on both the blue and red side of the systemic emission over a period of a 
few years. 

Conversely, if amplified fluctuations do not follow Keplerian motion, then it is re- 
markable that they are all located on the approaching side of the disk. I.e., individual 
maser features on the blue side are brighter and wider, > 50% in flux and about a factor 
of two in line width (Hashick et al. 1994), than their red counterparts. Also, fluctua- 
tions driven by turbulent motions have typical wavelengths of the order of H = 10 15 cm 
(Stone et al. 1996), while the height of the disk is generally much smaller than the extent 
Ar = 8 x 10 16 cm over which maser amplification takes place. So a large number number 
of fluctuations is expected to contribute to the systemic maser emission, on both the 
blue and red side. Hence, the asymmetry between red and blue systemic emission that 
is present in the Hashick et al. cumulative signal is unlikely to be random. In this work, 
these observations are explained by general relativity and saturation. Clearly, it is neces- 
sary to extend the Hashick et al. procedure for epochs upto the present. Note in this the 
respective integrals under the blue and red parts of the of Hashick et al. systemic maser 
emission. The surface area underneath the red systemic side is somewhat larger than 
that of the blue side. This implies that the total flux, i.e., bright features plus pedestal 
emission, on the blue systemic is set by a pumping rate and/or background radiation 
source that is not larger than that of the red side. Together with the much stronger and 
wider (saturated) maser features on the blue side this points to competitive gain among 
different maser beams on the blue systemic side (Alcock & Ross 1985). 

The strongest assumption that has been made in this work is that of constant excita- 
tion conditions. Still, models that invoke X-rays to heat the gas and drive ion- molecule 
reactions that lead to (warm) water favor oblique illumination, which yields constant 
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excitation condition when viewed edge-on (Neufeld et al. 1994). Hence, excitation condi- 
tions may indeed be fairly constant. Future work will concentrate on a detailed analysis 
of line profiles and excitation conditions in the presence of saturation. In addition, the 
impact of beam crossing (Spaans & van Langevelde 1992) on the shape of the line profile 
will be investigated then. 

Accretion disks around more massive black holes are desirable to further test general 
relativistic effects. Conversely, the effects mentioned here should be (statistically) smaller 
in lower mass systems like NGC 1068 (M ~ 10 7 Af Q ). Also, water maser galaxies much 
farther away than NGC 4258 do not allow the motions of individual maser spots to be 
followed. However, the observation of asymmetries in the systemic maser emission may 
then still lead to a (model dependent) determination of the black hole mass. Finally, for 
NGC 4258, gravitational redshift would not lead to observational effects if Ar < 0.08i?o- 
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